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between  the  two  ground  states.  In  most  cases 
(Fig.  4D),  one  of  the  modes  clearly  dominates. 
However,  in  some  cases  (below  10%),  both  con¬ 
figurations  are  almost  equally  populated  at  the 
same  time.  The  simultaneous  observation  of 
two  modes  might  be  due  to  excitations  like  the 
formation  of  spin  domains.  Hence,  for  the  first 
time,  our  system  allows  detailed  studies  on  the 
mechanism  of  symmetry-breaking  in  large-scale 
“magnetic”  systems. 

A  different  aspect  of  our  system  is  revealed 
when  one  recalls  that  the  spin  configurations  ac¬ 
tually  correspond  to  local  phases  of  a  BEC  at  dif¬ 
ferent  sites  of  the  triangular  lattice.  This  provides 
new  insight  into  unconventional  superfluidity 
(21,  22):  For  all  phases  except  the  ferromagnetic 
one,  the  state  corresponds  to  a  superfluid  at  non¬ 
zero  quasi-momentum  and,  thus,  nontrivial  long- 
range  phase  order.  Moreover,  the  observed  spiral 
configurations  spontaneously  break  time-reversal 
symmetry  by  showing  circular  bosonic  currents 
around  the  triangular  plaquettes  of  the  lattice. 
Clockwise  and  counterclockwise  currents  are 
found  in  a  staggered  fashion  from  plaquette  to 
plaquette.  These  resemble  the  currents  of  the  stag¬ 
gered  flux  state  thought  to  play  a  role  in  explain¬ 
ing  the  pseudogap  phase  of  high-temperature 
cuprate  superconductors  (23). 

Our  results  demonstrate  the  realization  of  a 
simulator  for  classical  magnetism  in  a  triangular 


lattice.  Let  us  point  out  here  that  these  results 
are  obtained  simply  by  using  spinless  bosons. 
Due  to  the  high  degree  of  controllability,  we  suc¬ 
ceeded  in  observing  all  of  the  various  magnetic 
phases  and  phase  transitions  of  first  and  second 
order,  as  well  as  frustration-induced  spontaneous 
symmetry-breaking.  It  now  becomes  possible  to 
quench  systems  on  variable  time  scales  from  fer¬ 
romagnetic  to  antiferromagnetic  couplings  and  to 
study  the  complex  relaxation  dynamics.  Further¬ 
more,  extending  the  studies  to  the  strongly  corre¬ 
lated  regime  promises  to  give  a  deeper  insight 
into  the  understanding  of  quantum  spin  models 
and  spin-liquid-like  phases  (24). 
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Visualizing  Individual  Nitrogen 
Dopants  in  Monolayer  Graphene 

Liuyan  Zhao/  Rui  He/  Kwang  Taeg  Rim/  Theanne  Schiros/  Keun  Soo  Kim/  4  Hui  Zhou/ 
Christopher  Gutierrez/  S.  P.  Chockalingam/  Carlos  ].  Arguello/  Lucia  Palova/ 

Dermis  Nordlund/  Mark  S.  Hybertsen/  David  R.  Reichman/  Tony  F.  Heinz/  7  Philip  Kim/ 
Aron  Pinczuk/  8  George  W.  Flynn/  Abhay  N.  Pasupathy1* 

In  monolayer  graphene,  substitutional  doping  during  growth  can  be  used  to  alter  its  electronic 
properties.  We  used  scanning  tunneling  microscopy,  Raman  spectroscopy,  x-ray  spectroscopy, 
and  first  principles  calculations  to  characterize  individual  nitrogen  dopants  in  monolayer 
graphene  grown  on  a  copper  substrate.  Individual  nitrogen  atoms  were  incorporated  as  graphitic 
dopants,  and  a  fraction  of  the  extra  electron  on  each  nitrogen  atom  was  delocalized  into  the 
graphene  lattice.  The  electronic  structure  of  nitrogen-doped  graphene  was  strongly  modified 
only  within  a  few  lattice  spacings  of  the  site  of  the  nitrogen  dopant.  These  findings  show  that 
chemical  doping  is  a  promising  route  to  achieving  high-quality  graphene  films  with  a  large 
carrier  concentration. 


Substitutional  doping  is  a  powerful  way  to 
tailor  material  properties  at  the  nanoscale 
(1),  and  it  might  be  expected  to  have  fun¬ 
damentally  different  consequences  when  used  to 
alter  the  electronic  properties  of  inherently  two- 
dimensional  (2D)  materials  such  as  graphene  (2). 
In  recent  years,  several  experimental  techniques 
have  been  developed  to  dope  the  carbon  lattice 
(3-9).  These  include  methods  applied  during 
growth  of  large-area  graphene  films  (4,  5,  7,  8) 
and  ways  to  modify  the  material  after  growth 
(3,  6),  as  well  as  a  new  one-pot  procedure  to 
produce  highly  doped,  few-layer  graphitic  struc¬ 


tures  (9).  Several  characterization  techniques 
including  x-ray  photoemission  (10, 11),  Raman 
spectroscopy  (3,  5,  7)  and  transmission  electron 
microscopy  (12)  have  been  used  to  analyze  the 
effect  of  the  doping  process  in  graphene.  How¬ 
ever,  a  microscopic  understanding  of  the  atomic 
and  low-energy  electronic  structure  induced  by 
the  substitutional  doping  process  in  monolayer 
graphene  is  lacking.  In  our  experiment,  we  used 
the  atomic-resolution  imaging  capabilities  of 
the  scanning  tunneling  microscope  to  understand 
the  local  structure  in  the  vicinity  of  a  N  dopant  in 
monolayer  graphene  and  spectroscopic  imaging 


to  measure  the  density  of  states  and  carrier  concen¬ 
tration  at  the  nanoscale.  We  also  report  scanning 
Raman  spectroscopy,  core-level  x-ray  spectroscopy 
and  first  principles  calculations  to  characterize 
the  effect  of  N-doping  on  the  graphene  film. 

We  grew  N-doped  graphene  films  using 
chemical  vapor  deposition  on  copper  foil  sub¬ 
strates  (13)  in  a  quartz  tube  furnace  [sample 
growth  details  in  supporting  online  material 
(SOM)  I  (a)].  Each  graphene-covered  foil  was 
divided  into  three  parts  and  immediately  studied 
using  Raman  spectroscopy,  x-ray  spectroscopy, 
and  scanning  tunneling  microscopy  (STM)  mea¬ 
surements.  Samples  have  been  characterized  as 
grown  on  copper  foils  as  well  as  after  transfer  to 
an  insulating  Si02  substrate  [sample  transfer  and 
preparation  details  in  SOM  I,  (b)  to  (d)]. 

Raman  spectra  (Fig.  1A)  were  taken  over 
large  areas  (20-pm  spot  size)  of  the  N-doped 
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samples  on  Si02  substrates  and  on  a  pristine 
graphene  sample  for  reference.  Raman  spectra 
taken  on  copper-foil  substrates  showed  identical 
features,  along  with  a  luminescence  background 
from  the  copper  [see  fig.  SI  and  SOM  II  (a)].  The 
pristine  graphene  spectrum  showed  the  G  and 
2D  resonances  at  1585.5  and  2684.5  cm-1  (14); 
little  or  no  signal  at  the  D  band  at  1348.5  cm-1 
was  observed,  indicating  a  low  defect  density.  In 
contrast,  the  doped  samples  exhibit  strong  D  and 
D'  resonances  ( 6 , 15-17),  as  well  as  changes  in 
the  intensity  and  position  of  the  G  and  2D  bands 


(18, 19).  The  Raman  spectra  on  N-doped  graphene 
were  consistent  with  the  presence  of  dopants  or 
defects  that  modify  the  free  carrier  concentration 
while  preserving  the  basic  structural  properties  of 
the  graphene  sheet. 

We  confirmed  the  presence  of  N  dopants  in 
the  graphene  lattice  by  performing  x-ray  photo¬ 
electron  spectroscopy  (XPS)  and  near-edge  x-ray 
absorption  fine  structure  (NEXAFS)  spectrosco¬ 
py  at  the  nitrogen  K-edge  for  a  pristine  and  NG1 0 
doped  sample  (Fig.  IB)  [see  SOM  I  (e)].  The 
addition  of  NH3  in  the  graphene  growth  re¬ 


sulted  in  sharp  peaks  at  400.7  and  408  eV  in  the 
NEXAFS  spectrum,  corresponding  to  ls-to-7i* 
and  -c*  transitions,  respectively,  for  a  single  mo¬ 
lecular  species.  The  sharpness  and  strong  polar¬ 
ization  dependence  of  the  peaks  indicates  that  this 
species  has  well-defined,  in-plane  N-C  bonds  in 
the  graphene  lattice.  Based  on  previous  studies 
(20-22),  this  peak  can  be  assigned  to  sp2-bonded 
graphitic  nitrogen  with  three  C  neighbors.  Rep¬ 
resentative  Nls  XPS  shows  a  higher-binding 
energy  component  (black  line)  for  the  NH3- 
grown  samples,  indicating  the  formation  of  more 


Fig.  1.  Raman,  XPS,  and  STM  spec¬ 
tra  of  N-doped  graphene.  (A)  Raman 
spectra  taken  at  pristine  and  N-doped 
(four  different  doping  levels)  graphene 
on  a  Si02/Si  substrate  (20-pm  laser 
spot)  showing  systematic  changes  in 
the  spectra  as  a  function  of  doping. 

A.U.,  arbitrary  units.  (B)  NEXAFS  (to¬ 
tal  electron  yield)  of  pristine  and  N- 
doped  graphene  on  copper  foil  at  the 
N  K-edge.  N-doping  results  in  a  new 
peak  in  the  spectrum  at  400.7  eV  due 
to  graphitic  N.  XAS,  x-ray  absorption  spectroscopy.  (Inset)  XPS  data  for  pristine 
and  NG10  graphene,  showing  a  higher-binding  energy  component  (black 
arrow)  for  the  doped  sample.  (0  Large-area  STM  image  of  N-doped  graphene 


-  2.0A 


1400  1600 

Raman  shift  (cm-1 


2600  2800 


395  400  405  410  415  420 

Photon  energy  (eV) 


2.0  A 


on  copper  foil  showing  the  presence  of  numerous  pointlike  dopants  (such  as 
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0/bias  =  0.8  V,  /set  =  0.8  nA). 


A 


C 


1- 0.34% 

2-  0.35% 

3-  0.34% 

4-  0.30% 

5-  0.30% 

6-  0.27% 

7-  0.27% 

8-  0.23% 


r  -  Distance  from  dopant  (nm) 


Fig.  2.  STM  imaging  of  nitrogen  dopants.  (A)  STM  image  of  the  most  common 
doping  form  observed  on  N-doped  graphene  on  copper  foil,  corresponding  to  a 
single  graphitic  N  dopant.  (Inset)  Line  profile  across  the  dopant  shows  atomic 
corrugation  and  apparent  height  of  the  dopant  (l/bias  =  0.8  V,  /set  =  0.8  nA).  (B) 
Simulated  STM  image  of  graphitic  N  dopant  (l/bias  =  0.5  V),  based  on  DFT 
calculations.  Also  superposed  is  a  ball-and-stick  model  of  the  graphene  lattice 
with  a  single  N  impurity.  (C)  STM  image  of  N-doped  graphene  on  copper  foil 
showing  14  graphitic  dopants  and  strong  intervalley  scattering  tails.  (Inset)  FFT 
of  topography  shows  atomic  peaks  (outer  hexagon)  and  intervalley  scattering 


peaks  (inner  hexagon,  indicated  by  red  arrow)  (l/bias  =  0.8  V,  /set  =  0.8  nA).  (D) 
Spatial  distribution  of  N-N  distances  from  eight  samples  on  copper  foils  with 
different  N  concentrations.  The  distributions  are  all  fit  well  by  a  quadratic  power 
law  (expected  error  bands  in  gray)  over  all  length  scales  indicating  that  N 
dopants  incorporate  randomly  into  the  graphene  lattice.  (E)  dl/dV  curves  taken 
on  a  N  atom  (bottom)  and  on  the  bright  topographic  features  near  the  nitrogen 
atom  on  N-doped  graphene  on  copper,  offset  vertically  for  clarity.  The  top  curve 
is  the  dl/dV  spectrum  taken  ~2  nm  away  from  the  dopant.  (Inset)  Positions 
where  the  spectra  were  taken  (l/bias  =  0.8  V,  /set  =  1.0  nA). 
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electronegative  N-C  bonds  in  the  graphene  lattice 
as  opposed  to  at  edges  or  defects  (gray  line).  This 
higher-binding  energy  peak  is  generally  consid¬ 
ered  a  signature  of  graphitic  N  in  studies  of  mod¬ 
ified  carbon  films  (6,  21,  23,  24). 

For  STM  measurements,  the  graphene-coated 
substrates  were  transferred  in  ambient  atmosphere 
soon  after  preparation  into  an  ultra-high-vacuum 
(UHV),  low-temperature  instrument  that  is  capa¬ 
ble  of  picometer  resolution.  The  treated  copper 
substrates  were  degassed  in  UHV  at  a  temper¬ 
ature  of  350°C  for  several  hours  before  STM  ex¬ 
periments  were  performed.  Although  many  areas 
of  the  copper  foil  displayed  a  rough  topography 
because  of  polycrystallinity,  we  occasionally  found 
large  areas  where  atomically  flat  terraces  can  be 
observed.  The  NG10  samples  were  optimal  for 
STM  studies,  because  there  is  a  reasonable  prob¬ 
ability  of  observing  such  flat  and  clean  areas 
along  with  the  dopant  features  described  below; 
our  detailed  STM  measurements  described  below 
are  performed  on  several  such  samples  between 
50  and  77  K. 


A  representative  STM  topography  of  N-doped 
graphene  on  a  large  copper  terrace  (Fig.  1C)  most¬ 
ly  revealed  the  pristine  graphene  lattice,  but  the 
image  also  showed  several  bright  objects  (nearly 
identical  in  appearance)  with  lateral  dimensions 
of  a  few  atomic  spacings,  consistent  with  previous 
measurements  on  doped  films  (9,  25).  These 
bright  objects  were  not  seen  in  pristine  graphene 
films  (see  below),  leading  us  to  associate  them 
with  N  dopants.  By  simply  counting  up  the  num¬ 
ber  of  bright  features  and  associating  each  with 
a  single  N  dopant,  we  arrived  at  a  N  doping  con¬ 
centration  per  carbon  atom  of  0.34%.  Measure¬ 
ments  across  several  similar  samples  yielded  N 
concentrations  between  0.23  and  0.35%. 

A  close-up  topography  of  one  of  these  dop¬ 
ing  features  (Fig.  2 A)  revealed  three  bright  spots 
forming  a  triangle.  The  distance  between  the  bright 
spots  was  equal  to  the  graphene  lattice  constant 
(2.5  A).  Far  away  from  the  doping  feature,  the 
honeycomb  lattice  of  graphene  was  recovered. 
Following  the  honeycomb  lattice  through  the  dop¬ 
ing  structure  situated  each  of  the  three  bright 


spots  on  a  C  atom  of  the  same  sublattice.  An 
STM  line  scan  through  the  dopant  (Fig.  2 A, 
inset)  yielded  a  maximum  apparent  out-of-plane 
height  of  0.6  ±  0.2  A  and  would  be  consistent 
with  the  N  atom  substituting  for  a  C  atom  in  the 
plane  of  the  graphene.  From  sampling  areas  on 
several  different  N-doped  graphene  samples,  we 
observe  that  more  than  90%  of  the  dopants  are 
of  this  form.  The  observed  STM  image  closely 
matched  our  simulated  STM  image  (Fig.  2B) 
computed  from  the  local  density  of  states  for 
graphitic  doping  where  one  N  atom  replaces  a 
single  C  atom,  as  well  as  other  recent  calculations 
(26).  Our  STM  image  simulations  were  carried 
out  using  the  Tersoff-Hamann  approach  (27)  based 
on  density  functional  theory  (DFT)  calculations 
in  the  local  density  approximation  for  exchange 
and  correlation,  with  a  plane-wave  basis  set  and 
norm-conserving  pseudopotentials  as  implemented 
in  the  Quantum-ESPRESSO  package  (28)  [see 
SOM  II  (e)].  The  calculated  STM  image  shows 
the  brightest  features  (red  spots)  on  the  nearest 
neighbor  C  atoms  to  N  (shown  in  blue).  Visible 
features  extended  to  several  lattice  spacings  from 
the  N,  and  the  overall  triangular  symmetry  of  the 
image  matched  the  experiment. 

The  STM  observations  show  that  the  majority 
of  the  doping  occurs  via  graphitic  substitution  in 
these  samples  [see  SOM  II  (b)  and  fig.  S2  for 
examples  of  other  doping  forms  seen].  Larger- 
area  STM  images  of  the  samples  (such  as  Fig. 
2C)  located  14  dopants,  all  in  the  graphitic  form, 
which  is  consistent  with  the  NEXAFS  results. 
Apart  from  the  local  structure  around  each  dop¬ 
ant,  we  also  see  long  “tails”  around  each  dopant. 
These  arise  from  intervalley  electron  scattering 
induced  by  the  N  dopant,  and  similar  features  have 
been  seen  before  in  other  experiments  (29,  30). 
The  fast  Fourier  transform  (FFT)  of  Fig.  2C  (inset) 
shows  evidence  for  this  strong  scattering.  The 
FFT  shows  two  sets  of  points  arranged  in  hexa¬ 
gons.  Whereas  the  outer  hexagon  of  points  in  the 
FFT  corresponds  to  the  atomic  lattice,  the  inner 
hexagon  in  this  FFT  arises  from  the  strong  inter¬ 
valley  scattering  induced  by  N  dopants  (29,  30). 

Large-area  STM  images,  such  as  those  in  Figs. 
1C  and  2C,  can  give  us  information  on  whether 
clustering  of  dopants  occurs  during  growth.  To 
analyze  the  possibility  of  dopant  clustering,  we 
used  STM  images  to  calculate  all  of  the  N-N 
distances  for  each  sample.  Figure  2D  shows  a 
histogram  of  these  distances  for  several  different 
samples  on  a  logarithmic  scale.  A  true  random 
distribution  of  dopants  would  imply  a  quadratic 
distribution  of  dopant-to-dopant  distances,  which 
we  observed  in  all  of  our  samples  down  to  length 
scales  of  a  few  lattice  constants.  Nearby  dopants 
preferred  to  incorporate  on  the  same  sublattice  of 
graphene.  Indeed,  all  of  the  dopants  shown  in 
Fig.  2C  are  located  on  the  same  sublattice,  and 
evidence  for  this  phenomenon  is  seen  in  other 
samples  as  well  [see  SOM  II  (c)  and  fig.  S3]. 

We  performed  detailed  STM  spectroscopy 
measurements  of  the  differential  conductance 
dl/dV  (Fig.  2E)  to  understand  the  effect  of  the 
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Fig.  3.  Spectroscopy  and  doping  in  N-doped  graphene.  (A)  Histogram  of  the  spatial  variation  of  the  Dirac 
point  for  N-doped  graphene  on  copper  foil  over  an  area  of  40  nm  by  40  nm  (l/bias  =  0.6  V,  /set  =  1.0  nA). 
(Inset)  Spatially  averaged  dl/dV spectrums  (black  line)  and  variation  (gray  band)  over  the  area.  The  orange 
arrow  indicates  the  position  of  the  Dirac  point.  (B)  Spatially  averaged  graphene  charge-carrier 
concentration  as  a  function  of  average  nitrogen  concentration  level  for  five  different  samples 
measured  on  N-doped  graphene  on  copper  foil.  (Inset)  Free  charge  carriers  per  nitrogen  atom  for 
each  of  these  samples.  (C)  Calculated  projected  density  of  states  near  the  Fermi  energy  for  a  1% 
doping  level.  The  gray  dashed  line  denotes  the  position  of  the  Fermi  energy.  (Inset)  Color  plot  of  the 
calculated  charge  density  around  the  N  dopant.  (D)  Averaged  dl/dV  spectrum  for  N-doped  graphene 
transferred  to  a  Si02/Si  substrate  (black  line)  and  spatial  variation  (gray  band).  The  Dirac  point  is  at  -330  meV 
(red  arrow)  (l/bias  =  1.0  V,  /set  =  0.1  nA).  (Inset)  Topography  shown  in  derivative  mode  to  enhance  atomic 
contrast. 


www.sciencemag.org  SCIENCE  VOL  333  1 9  AUGUST  201 1 


1001 


Downloaded  from  www.sciencemag.org  on  September  1 , 201 1 


REPORTS 

N  atoms  on  the  low-energy  electronic  structure 
of  the  graphene  film  (here,  dlldV  is  the  deri¬ 
vative  of  the  current  with  respect  to  the  voltage 
obtained  using  a  lock-in  amplifier).  The  most 
prominent  features  seen  in  all  of  the  curves, 
which  have  been  seen  previously  on  monolayer 
graphene  (31),  are  two  depressions  near  zero 
bias  and  near  -300  meV  relative  to  the  tip.  The 
zero-bias  feature  occurs  when  energies  are  too 
low  to  excite  the  optical  phonon  mode  in  graphene 
monolayers  that  can  enhance  the  tunneling  cur¬ 
rent  (31).  The  -300  meV  feature  occurs  near 
the  Dirac  point  where  the  density  of  states  of 
graphene  is  low  (31-33).  We  thus  associated  the 
dip  seen  at  -270  mV  in  Fig.  2E  with  the  Dirac 
point  for  an  electron-doped  graphene  layer.  Mea¬ 
surements  of  the  Dirac  point  over  a  40-  by  40-nm 
area  of  an  NG10  sample  show  a  distribution  of 
values  (Fig.  3A).  We  used  the  energy  of  the 
Dirac  point  measured  at  each  position  of  a  sample 
to  convert  to  a  charge-carrier  density  at  that  loca¬ 
tion  using  the  relation  for  the2ideal  graphene  band 
structure — namely,  n  =  ^ .  2  (where  n  is  the 
charge-carrier  density,  ED  is  the  Dirac  point 
energy,  and  h  is  Planck’s  constant  h  divided  by 
2ri).  Taking  a  value  for  the  Fermi  velocity  vF  = 
106  m/s,  we  arrive  at  an  average  charge-carrier 
density  of  (5.42  ±  0.83)  x  1012  electrons  per  cm2. 
We  can  compare  this  to  the  N-doping  concentra¬ 
tion  calculated  by  simply  counting  the  number  of 
observed  dopants  in  the  area  where  the  measure¬ 
ments  are  performed.  The  observed  N  doping  in 
this  area  corresponds  to  0.34%  N  atoms  per  C 
atom  or,  equivalently,  a  N  atom  density  of  1.3  x 
1013cm-2.  Together  with  the  charge-carrier  density 
measured  by  STM  spectroscopy,  this  result 
implies  that  each  graphitic  N  dopant  contributes 
(on  average)  -0.42  ±  0.07  mobile  carriers  to  the 
graphene  lattice.  We  have  performed  such  de¬ 
tailed  spectroscopic  measurements  across  differ¬ 
ent  samples  with  N  atom  concentrations  varying 
from  0.23  to  0.35%  (Fig.  3B)  and  found  a  strong 
correlation  between  the  number  of  N  dopants 
and  the  extracted  charge-carrier  density.  An  in¬ 
dependent  measurement  of  the  carrier  concen¬ 
tration  based  on  the  shift  in  the  position  of  the  G 
peak  in  the  Raman  spectrum  induced  by  doping 
(18, 19)  was  5  ±  1 .5  x  1 012  cnT2  [see  SOM  II  (d) 
and  fig.  S4]. 

Our  DFT  calculations  (28)  provide  insight  to 
the  electrostatic  balance  between  the  N  dopants 
and  free  carriers  in  the  graphene  sheet.  Focusing 
on  a  single,  graphitic  N  dopant,  the  projected  N 
density  of  states  (pDOS)  on  the  7i-system  (Fig. 
3C)  revealed  a  resonance  caused  by  the  N  pz  or¬ 
bital  centered  0.3  eV  above  the  Fermi  energy  (EF). 
The  pDOS  for  the  C  nearest  neighbor  exhibited  a 
shoulder  caused  by  its  electronic  coupling  to  the 
N,  with  a  reduced  shoulder  on  the  next  nearest 
neighbor.  The  occupied  fraction  of  these  reso¬ 
nances  represented  the  localized  charge  near  the 
N  centers.  The  balancing  charge  went  to  the  rest 
of  the  71  states.  The  Dirac  point  still  appeared  in 
the  pDOS,  shifted  to  below  E¥,  as  also  seen  in 
our  total  DOS  and  in  other  recent  calculations 


(26,  34).  We  studied  N  concentrations  (from  0.6 
to  5.6%)  by  varying  the  super  cell  size  for  a 
single  dopant.  Using  the  Dirac  point  shift  in  the 
total  DOS  as  a  measure  of  free  electron  concen¬ 
tration,  we  estimate  (using  the  formula  above) 
that  50  to  70%  of  electrons  are  delocalized,  with 
oscillations  indicative  of  electronic  interference 
effects  in  the  simulation  cells  that  we  used  [see 
SOM  II  (e)  and  fig.  S5].  A  value  of  -  60%  would 
be  consistent  with  the  experiment  values  shown 
in  Fig.  3B. 

To  draw  comparisons  with  other  measure¬ 
ments  of  doping  in  graphene  films,  it  is  important 
to  understand  the  effect  of  the  copper-foil  sub¬ 
strate  on  the  carrier  concentration  via  charge  trans¬ 
fer  (35,  36),  as  well  as  by  changing  the  charge 
screening  in  the  graphene  layer.  We  performed 
STM  measurements  of  pristine  graphene  films  on 
copper  foils  and  estimated  that  the  doping  induced 
by  the  copper-foil  substrate  into  the  graphene 
film  is  <1012  electrons/cm2  [see  SOM  II  (f)  and 
fig.  S7].  The  copper  substrate  can  also  modify 
the  charge  screening  length  in  the  graphene  film. 
We  studied  doped  graphene  films  transferred  to  a 
Si02  dielectric  substrate  (37),  a  process  that  left 
residue  on  or  below  the  graphene  surface.  Typical 
STM  and  atomic  force  microscopy  images  dis¬ 
played  surface  roughness  of  a  few  nanometers, 


but  occasionally  we  found  small  regions  of  the 
sample  where  the  graphene  honeycomb  lattice 
could  be  resolved.  Figure  3C  shows  the  average 
dl/dV  spectrum  taken  over  one  such  clean  area 
(the  STM  image  is  shown  in  the  inset  in  the  de¬ 
rivative  mode  to  remove  the  overall  roughness 
of  the  substrate  and  enhance  the  atomic  contrast). 
The  overall  features  in  the  spectrum  are  preserved, 
but  now  Ed  =  -330  ±  20  meV  (n  =  8  x  1 012cnf 2). 
The  transfer  processes,  as  well  as  the  Si02  sub¬ 
strate  itself,  introduce  an  unknown  doping  con¬ 
centration  into  the  film;  thus,  it  is  not  possible  to 
directly  compare  the  doping  level  of  the  two  sam¬ 
ples.  However,  the  spectrum  of  doped  graphene 
on  Si02  is  broadly  consistent  with  the  spectrum 
of  doped  graphene  on  copper  foil,  indicating  the 
absence  of  a  strong  hybridization  between  the 
graphene  and  underlying  copper. 

The  local  electronic  structure  around  a  N  dop¬ 
ant  also  affects  the  electronic  nature  of  the  doped 
graphene  film.  In  Fig.  4A,  we  show  spectra  ob¬ 
tained  on  and  far  away  (-2  nm)  from  a  dopant 
atom.  Although  the  overall  features  of  the  spectra 
were  preserved  on  the  N  atom,  the  electron-hole 
asymmetry  in  the  local  density  of  states  was 
much  stronger  on  the  N  atom,  in  accord  with  DFT 
calculations  (Fig.  3B).  The  amount  of  scattering 
introduced  by  the  dopant  can  be  assessed  by 


Fig.  4.  Spectroscopic  mapping  around  a  single  N  dopant  (A)  dl/dV  spectra  taken  on  the  N  atom  (solid  line) 
and  far  away  from  the  N  atom  (dashed  line)  for  N-doped  graphene  on  copper  foil  (l/bias  =  0.8  V,  /set  =  1.0  nA). 
(B)  STM  spectroscopic  maps  taken  in  the  vicinity  of  a  single  N  dopant  for  N-doped  graphene  on  copper  foil 
(dopant  position  indicated  by  a  dot)  at  different  energies  with  N  positions  marked  by  red  diamonds  (l/bias  = 
0.8  V,  /set  =  1.0  nA).  (C)  Radially  averaged  differential  conductance  as  a  function  of  the  distance  from  the 
N  atom  at  different  energies,  normalized  to  unity  at  distances  far  away  from  the  N  atom.  The  fits  are  to  an 
exponential  decay  function.  (D)  Extracted  decay  lengths  for  different  energies  (squares)  and  the  ratio  of 
the  differential  conductance  on  the  N  site  to  the  background  (triangles)  for  different  energies.  The  decay 
lengths  are  ~7  A  for  all  energies. 
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STM  imaging  (29,  32,  38).  We  performed  STM 
spectroscopic  mapping  in  a  2.5-nm  area  around 
one  graphitic  N  dopant  (STM  junction  param¬ 
eters  are  V=  0.8  V,  1=  1  nA).  Shown  in  Fig.  4B 
are  a  subset  of  these  maps,  acquired  at  bias  volt¬ 
ages  (applied  to  the  sample  relative  to  the  tip) 
from  -0.78  to  0.78  V.  The  maps  did  not  show 
much  contrast  at  high  positive  bias,  but  the  local 
DOS  around  the  N  atom  was  strongly  suppressed 
at  energies  below  the  EF.  The  local  DOS  recov¬ 
ered  its  background  value  within  a  few  lattice 
constants  of  the  dopant  atom.  We  plotted  the 
radial  distribution  of  the  dl/d  V  intensity  (Fig.  4C) 
from  the  maps  in  Fig.  4B  as  a  function  of  distance 
from  the  N  atom,  normalizing  the  background 
value  of  the  dl/dV  to  unity  for  each  energy.  The 
intensity  of  the  spectral  weight  changes  caused  by 
the  N  dopant  were  energy-dependent,  but  the  de¬ 
cay  lengths  were  ~7  A  for  all  energies  (Fig.  4D). 
This  result  indicates  that  the  electronic  perturba¬ 
tion  caused  by  a  nitrogen  dopant  is  localized  near 
the  dopant  atom,  which  is  confirmed  in  large-area 
dl/dV  maps,  and  seen  in  the  calculated  charge 
distribution  [see  SOM  II  (e)  and  fig.  S5]  and 
simulated  STM  image  in  Fig.  2B  (26). 
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Guided  Growth  of  Millimeter-Long 
Horizontal  Nanowires  with 
Controlled  Orientations 

David  Tsivion,1  Mark  Schvartzman,1  Ronit  Popovitz-Biro,2  Palle  von  Huth,2  Ernesto  Joselevich1* * 

The  large-scale  assembly  of  nanowires  with  controlled  orientation  on  surfaces  remains  one 
challenge  preventing  their  integration  into  practical  devices.  We  report  the  vapor-liquid-solid 
growth  of  aligned,  millimeter-long,  horizontal  GaN  nanowires  with  controlled  crystallographic 
orientations  on  different  planes  of  sapphire.  The  growth  directions,  crystallographic  orientation, 
and  faceting  of  the  nanowires  vary  with  each  surface  orientation,  as  determined  by  their 
epitaxial  relationship  with  the  substrate,  as  well  as  by  a  graphoepitaxial  effect  that  guides  their 
growth  along  surface  steps  and  grooves.  Despite  their  interaction  with  the  surface,  these 
horizontally  grown  nanowires  display  few  structural  defects,  exhibiting  optical  and  electronic 
properties  comparable  to  those  of  vertically  grown  nanowires.  This  paves  the  way  to  highly 
controlled  nanowire  structures  with  potential  applications  not  available  by  other  means. 


Semiconductors  with  controlled  structures  are 
at  the  core  of  the  most  advanced  technolo¬ 
gies  in  our  everyday  life,  but  defects  induced 


by  stress  during  growth  often  degrade  the  elec¬ 
tronic  and  optical  properties.  For  example,  single¬ 
crystal  gallium  nitride  was  first  produced  in  1969 


by  epitaxial  growth  on  sapphire  (1),  but  it  took 
three  decades  of  research  to  reduce  its  defect  con¬ 
centration  before  it  became  the  basis  of  the 
blue  light-emitting  diodes  (LEDs)  and  violet  lasers 
that  enabled  outdoor  TV  screens  and  the  Blu- 
ray  disk.  The  vapor-liquid-solid  (VLS)  growth 
method,  first  described  in  1964  (2),  has  attracted 
overwhelming  renewed  attention  during  the  past 
decade  (3),  as  a  means  of  producing  stress-free 
single-crystal  semiconductor  nanowires  with  un¬ 
paralleled  electronic  and  optical  properties  suitable 
forultraminiaturized  electronics  (4),  optoelectronics 
(5),  renewable  energy  (6),  and  chemical  and  bio¬ 
logical  sensing  (7).  Yet  the  large-scale  assembly 
of  horizontal  nanowires  with  controlled  orienta¬ 
tion  on  surfaces  remains  a  challenge  toward  their 
integration  into  practical  devices.  To  this  end,  sev¬ 
eral  assembly  processes  have  been  devised,  includ¬ 
ing  the  use  of  liquid  flows  (8),  electric  fields  (9), 
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